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Abstract 

We use the Wilson-Fisher e expansion to study quantum critical behavior in gauged 
Yukawa matrix field theories with weak quenched disorder. We find that the resulting 
quantum critical behavior is in the universality class of the pure system. As in the 
pure system, the phase transition is typically first order, except for a limited range 
of parameters where it can be second order with computable critical exponents. Our 
results apply to the study of two-dimensional quantum antiferromagnets with weak 
quenched disorder and provide an example for fluctuation-induced first order phase 
transitions in circumstances where naively none is expected. 



One of the most remarkable features of quantum spin systems is their relationship to 
gauge theories. This connection, which was originally used to study chiral symmetry breaking 
in QCD in the strong coupling limit [J], remains a powerful tool. Recently this analogy 
has been used to prove that certain gauge theories break chiral symmetry in the strong 
coupling limit [2-4]. Also, it has been further suggested || that the underlying analogy 
can be used in a much broader sense to account for the quasi-particle spectrum and other 
infrared features of two-dimensional (2D) quantum antiferromagnets and three-dimensional 
non-Abelian gauge theories. Using the quantum chromodynamics (QCD) terminology, the 
quantum phase transition with a chiral symmetry breaking pattern is then characterized 
by the number of flavors (Np) and colors (Nc) of quarks in the gauge theory. In the 
antiferromagnet the rank of the spin algebra and the size of its representation play the same 
role as the number of flavors and colors respectively. 

The study of 2D quantum antiferromagnets, particulary in connection with high-Tc su- 
perconductivity, is a rapidly developing subject of great current interest |J. An important 
question that naturally arises is the stability of the antiferromagnetic long-range order (LRO) 
in the presence of quenched disorder (QD). In fact, since the Neel temperature of planar high- 
T c superconductors, such as La 2 Cu0 4 , is extremely sensitive to impurities and defects, it 
is probably not correct to ignore the impurity effects even in high-quality samples used in 
laboratories. The critical behavior of low- dimensional systems with quenched disorder and 
an 0(A r )-vector magnetic order parameter has been addressed by a number of authors 0, [8| 
by using a Landau- Ginzburg- Wilson (LGW) Hamiltonian and applying the Wilson-Fisher 
(WF) e expansion ||. In particular, Boyanovsky and Cardy [f| have carried out a two- loop 
(double) e expansion in systems in which the impurities are correlated in dimensions and 
randomly distributed in d — dimensions. They take the full structure of the theory into 
account which, as a result of anisotropies, leads to highly nonlocal interactions. 

In this paper we shall study the stability of LRO in 2D quantum antiferromagnets with 
QD by using the WF e expansion and renormalization group (RG) techniques. We shall 
approach this problem by examining the critical behavior of the related zero temperature 
gauged Yukawa matrix field theories in three (Euclidean) dimensions in the presence of 
quenched disorder. This requires that the anisotropies be properly taken into account, as is 
done (see Ref. |§, e.g.) for the 0(iV)-vector model in the presence of QD. However we shall 
go one step further and also examine the stability of the theory against quantum fluctuations 
which can arise through the Coleman- Weinberg mechanism flQ| . Our results enable us to 
argue that below a certain critical doping (impurity concentration) A crit , and in a narrow 
region in parameter space, 2D doped quantum antiferromagnets belong to the universality 
class of gauged Yukawa matrix field theories which describe the pure system. We check 
this result through an explicit calculation of the critical exponent v and find that the Harris 
criterion Unjl is indeed satisfied. However, we find that a careful examination of the RG flows 



in the space of couplings indicates that LRO can be destroyed due to quantum fluctuations. 
This means that there are cases where there is an infrared (IR) stable fixed point and the 
Harris criterion is satisfied, but the phase transition is fluctuation-induced first order rather 
than second order, contrary to what is naively expected. We shall begin with a brief review 
of some recent work on pure quantum antiferromagnets and gauged Yukawa matrix theories. 
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The pure 2D (generalized) quantum antiferromagnet is defined by the Hamiltonian 



N%-1 



tfspi* = « E E J A (x)J A (v) (i) 

<x,y> A=l 

with < x,y > nearest neighbor sites x and y on a square lattice and the "spin" operators 
J A (x) in an irreducible representation of the SU(Np) Lie algebra 

~J A (x),J B (y)} = if ABC 6(x,y)J c (x) (2) 

When the representation at each site is a rectangular Young Tableau with m rows and Nc 
columns, it is convenient to represent the spin operators by the fermion bilinears 

N c N F 

J A {x) = E E ^L(*)zftM*) (3) 

a=l a,b=l 

The fermions have the anticommutator, 

{VVkO), tplp(y)} = S ab 8 a p8(x, y) (4) 

Constraints which project out the irreducible representation of the spin algebra are 

N F 



1>laWMx) ~ ~ Vx (5) 



^ Q/ g(x) obeys the U(Nc) Lie algebra, commutes with the Hamiltonian and acts as the gen- 
erator of gauge transformations with gauge group U (Nc) ■ 



The critical behavior of the antiferromagnet was examined by Read and Sachdev JL2 
using semiclassical methods. The only free parameters are the integers Nq and Np. Nc 
N F is the classical limit of large representations, where the classical Neel ground state is 
stable with the staggered spin order parameter 

N c 

^ab = {-l)^ Xl < E ^aa(x)lp ba (x) > (6) 
a=l 

On the other hand, the limit Np <C Nc is the quantum limit where fluctuations are important 
and the system is in a spin disordered state. For both Nc and Np large, they find a line of 
second order phase transitions in the (Nc, Np) plane at iV^ = const. • iVc where the constant 
is a number of order one. 

The relationship between the antiferromagnet and QCD is a very close one. There is an 
argument in ref. which maps the strong coupling limit of lattice QCD onto the quan- 
tum antiferromagnet. The lattice regularization of the QCD Hamiltonian uses staggered 
fermions [T3[. Since the order parameters are identical, the Neel ordered states of the an- 



tiferromagnet correspond to chiral symmetry breaking states of QCD. Thus, the infinite 
coupling limit of QCD is identical to the quantum antiferromagnet. A main difference be- 
tween QCD with finite coupling and the quantum antiferromagnet is that QCD allows a 
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fermion kinetic energy term while retaining gauge invariance, whereas in the antiferromag- 
net, the fermions are not allowed to move. One could, however, regard the corrections to the 
strong coupling limit of QCD as the addition of degrees of freedom and gauge invariant per- 
turbations in the quantum antiferromagnet which allow fermion propagation. In fact, it has 
even been suggested || that the additional degrees of freedom are generated dynamically. 

A common feature of both theories is that, aside from Np and Nq they have no free 
parameters. One could imagine adding operators of the sort that, if their coupling constants 
is varied, can induce the chiral phase transition. It is tempting to speculate that these 
transitions fall into a universality class which can take into account all such modifications, 
as long as they respect the symmetries of the theory. Restricting attention only to those 
operators which lead to a Lorentz invariant continuum limit, it was recently argued JT3J that 
the universality class is described by the 4 — e dimensional Euclidean field theory, 



S = I d 4 ~ e x 



itr(V0) • (V0) 



+ ^r(iv| (tr0) + Av tr0 ) + 4 tr ^ 

+^(7 ■ V + z// /2 ei7 • A + i// /2 e 2 7 ■ tiA 



Here the scalar field <fi is an Np x Np traceless Hermitean matrix. The 4-component spinor 
ip is an Np x Nq complex matrix and A^ is a U(Nq) gauge field. In four dimensions this 
model has Euclidean Lorentz invariance, C,P and T, discrete chiral symmetry, (ip — > 7 5, 0, 
cf) — > —(f)) and global SU(Np) flavor. (0) includes all operators which are marginal when 
D = 4. 



The evidence that ([?]) describes the universality class comes from previous work [15| 



where a similar model where gauge couplings are absent was studied. In Ref. jnj it was 
shown that the anomalous dimensions of operators computed in the model (|7[) with = 
were identical to leading order in 1/Nc and e to those of a 2 < D < 4 dimensional four- 
fermi theory. That a 4 — e dimensional Yukawa-Higgs theory has the same universal critical 
behavior as lower dimensional four-fermi theories with the same symmetries was originally 



suggested by Wilson ||16||. For the case Np = 1, where the chiral symmetry is discrete, 



higher order computations have been carried out |17|, fij|. These results, as well as those of 
lattice simulations, support the universality hypothesis |Tfl| . It was argued in [TJ| that (0) 
represents the universality class of lower-dimensional four-fermi theories with U(Nc) gauge 
invariance such that for a large range of values of (Nc, Np), the chiral phase transition is a 
fluctuation induced first order transition and when it is second order, critical exponents are 
in principle computable in the e expansion. 

To include the effects of QD amounts to adding new operators which can induce the chiral 
phase transition by varying the width, A, of the impurity-probability distribution. This is 
achieved by adding 

d d x{r + 5r(x))tr(f) 2 (x) (8) 
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to the aciton (|7|). Here 5r(x) is the impurity field and r is the mass of the field <f>{x) 
(renormalized to zero at the critical point). The quenched partition function can be written 
by using the replica trick ||19|| . Defining the impurity-probability distribution P[Sr] such that 
the replica average ((5r(x))) = and 

((6r(x)6r(x'))) = A5 d - td {x - x'), (9) 

a cumulant expansion yields the (replica) action for the gauged Yukawa matrix field theory 
with QD 

A 



Sqd = £ S « ~ % J ^ J ddx ' S<i ' ed ( x ~ x ') 



Etr^(x)Et4(x'). (10) 

a 13 

Here a, (3 = 1, . . . ,n are the replica indices and S a has the same form as S in ([?[) with a replica 
index a for each field operator and d — 4 — e. In ( fT0|) the impurities are correlated along 
€d dimensions, whereas in the remaining d — €d directions the impurity-correlation function 
(§) is zero for impurities at different points. It is important to note that the system is 
characterized by two dimensionalities, namely e and e^, which means that the perturbative 
expansion will, in general, involve a double e expansion. Such a double e expansion is 
necessary for a controlled perturbative description of the critical behavior (see Dorogovtsev 
in Refs. and §). 

To investigate the critical behavior we first compute the f3 functions for the model defined 
by flrop. To leading order, our calculations yield 

Nl + 7 2 2M-3 + 3 2 1 2 

A - -^ + ^ 1 2 + ^^ + ^ 2 2 + ^-^A (ii) 

P2 = - £9 * + W/ l9 * + ^ 9 >-8^/ + ^/ 92 > (12) 

a e 3 \ 2 3 2 N 2 F + 2N c N F -3 3 

^ = -2 2/ -16^^V^ ei2/ -8^ e22/ + 16iVciV| ^' (13) 

- £ c, ^-2^ 3 (14) 

- 2 ei 48?r2 e 1; (14j 

A. = + (15) 
(3 A = -~eA - 16A 2 + -^-yiA + -^^A, (16) 

where e = e + and a factor of 27r' i r((i/2), with o? = 4 — e, has been absorbed into the 
definition of A. Fixed points occur at the zeros of the f3 functions, /3(g*) = 0. The fixed 
points are IR stable if all the eigenvalues of the stability matrix, dj3i/dgj\ g=g *, are positive. 
Second order phase transitions are possible when the RG trajectories flow to the IR stable 
fixed point in the region allowed by general stability conditions. Otherwise, the only allowed 
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phase transition is a fluctuation induced first order one. Yamagishi |2(J formulated a criterion 



for this behavior. He showed that this occurs when the RG trajectory crosses the surface in 
coupling-constant space given by 

Vi(g,y,e) = 0, i = 1, 2, (17) 

where V l (g, y, e) = (4 - e) (g x + g 2 ) + f3 1 + (3 2 and V 2 = (4 - e) ( 5l + (N F /2)g 2 ) + & + (N F /2)(3 2 
depending on whether g 2 > or g 2 < 0, respectively. When the RG trajectories cross 
these surfaces two further conditions must be met. To ensure that the extremum is a local 
minimum, rather than a maximum, it is necessary that Di > for i — 1, 2 where D\ = 
(4-e)(/? 1 +/3 2 )+£ i A<9/%(/3 1 +/? 2 ) or D 2 = (4-e)(/5 1 + (iV F /2)/5 2 ) + S i /5^/%(A + (iV F /2)/5 2 ) 
depending on whether g 2 > or g 2 < 0. 

For this minimum to have lower free energy than the one at the trivial <fi = 0, it is 
necessary that the couplings at that scale obey 

gi + g 2 <0 or gi + N F g 2 /2<0. (18) 

By using Eqns. (11)-(16) we have numerically solved for the RG trajectories in the six- 
dimensional space of couplings. In Fig. 1 we have plotted the projection of some of the RG 
trajectories onto the g\ — g 2 plane for the case Nc = 2, N F = 13 and e = = 0.1. As in 
the pure case, some of the RG trajectories satisfy the Yamagishi conditions for fluctuation- 
induced first order behavior before reaching the IR stable fixed point. Two examples of this 
are the middle dashed and lower solid curves corresponding to cases with and without QD 
respectively . 

To guarantee the consistency of the above analysis based on the WF e expansion one 
needs to check for the Harris criterion [|l l| which requires the calculation of the critical 
exponent v given by 

"= (19) 

where 7 2 (g*) is the (mass) anomalous dimension. We have calculated 7 2 (g) in the presence 
of QD: 

Nl + 1 2Nl - 3 1 2 „ A 

72 = o 9i + — ^^92 + y - 4A 20 

' 6N F y 6N F y AN F V ; 

and checked for the conditition dv > 2 (the Harris criterion) at the fixed point (g*,g 2 )- 
Although we find that the Harris criterion is satisfied in the entire region < e, e < 1, 
an examination of the RG trajectories (see Fig.l) reveals that as in the pure case, except 
for a narrow range of parameters in the (Nc, N F ) plane, the phase transition is fluctuation 
induced first order. The occurence of this phenomenon depends on the parameters iV^ 
and N F (see Ref. [|rjj for details on a similar system without impurities) and typically 



there exists a critical impurity concentration, A crit ., above which the phase transition will 
be fluctuation induced first order. In the range of parameters where the phase transition is 
second order, the gauged Yukawa matrix theory with (weak) quenched disorder belongs to 
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the universality class of the pure system. Since the specific heat exponent remains neagative, 
this is expected and further supports the consistency of this work with previous results. The 
physical quantum spin j antiferromagnet corresponds to the parameter values Np = 2, 
Nc = 2j. With Nc = l(j = 1/2) the non-Abelian field disappears and the resulting theory 
has no IR fixed point, indicating a first order phase transition. The case Np = 2, Nc > 2 (i.e., 
j > 1) belongs to the asymptotically free regime and hence these antiferromagnets cannot 
be analyzed by our techniques. In this domain in the N C ,N F parameter space we speculate 
that confinement is associated with a nonperturbative IR fixed point of the gauge coupling. 
In this case these antiferromagnets are likely to exhibit second order phase transition. 

In conclusion, we have shown that long wavelength quantum fluctuations can drive first 
order phase transitions in the presence of weak quenched disorder when naively none is 
expected. It is important to note that we have studied the effects of QD on the critical 
behavior of the pure system by first analyzing the relevance or irrelevance of disorder at the 
critical fixed point of the pure system and then carrying out a thorough examination of the 
stability of the pure system against quantum fluctuations. On the other hand, one may also 
study the crossover from the pure fixed point to the disorder fixed point. This is an important 
question whose examination, however, is beyond the scope of the present paper. In general, 
detailed calculations [21] show that when the pure fixed point is unstable (i.e. vd < 2) the 



disorder fixed point is a critical fixed point with new critical exponents [In such instances the 
stability matrix typically has complex eigenvalues (see, e.g., Dorogovtsev in Ref. 0), which 
imply that the RG trajectories spiral into the fixed point basins of attraction, and give 
oscillatory corrections to scaling]. The results of this paper may help us better understand 
the nature of phase transitions and the underlying dynamics in low-dimensional systems 
with weak quenched disorder, such as those occuring in planar high-T^ superconductors. 

It is a pleasure to thank Michael Ma for invaluable discussions and correspondence 
which culminated in this work. I would also like to thank T.H. Hansson, E. Langmann 
and L.C.R. Wijewardhana for thoughtful comments and numerous discussions while this 
work was being completed. 

This work is supported in part by the United States Department of Energy under Grant 
No. DE-FG02-84ER40153 and by the Swedish Natural Science Research Council. 

References 

[1] J. Smit, Nucl. Phys. B175, 307 (1980); T. Banks et ai, Phys. Rev. D 15, 1111 (1977). 
[2] M. Salmhofer and E. Seiler, Commun. Math. Phys. 139, 395 (1991). 
[3] E. Langmann and G. Semenoff, Phys. Lett. 297B, 175 (1992). 

[4] M.C. Diamantini, P. Sodano, E. Langmann and G. Semenoff, Nucl. Phys. B406, 595 
(1993). 



[5] P. Beran, P. Poilblanc and R. B. Laughlin, preprint |cond-mat / 9505085 



7 



[6] V. Barzykin and D. Pines, Phys. Rev. B 52, 13585 (1995); S. Chakravarty, B. I. Halperin 
and D. R. Nelson, Phys. Rev. B 39, 2344 (1989). 



[7] A.B. Harris and T.C. Lubensky, Phys. Rev. Lett. 33, 1540 (1974); G. Grinstein and 
A. Luther, Phys. Rev. B 13, 1329 (1976); A. Aharony in Phase Transitions and Critical 
Phenomena, edited by C. Domb and M.S. Green (Academic Press, New York, 1976), Vol. 
VI, p. 357; T.C. Lubensky in Ill-Condensed Matter, edited by R. Balian, R. Maynard 
and G. Toulouse (North- Holland, Amsterdam, 1979); D.E. Khmelnitskii, Phys. Lett. 
67A, 59 (1978); S.N. Dorogovtsev, Phys. Lett. 76A, 169 (1980); T. R. Kirkpatrick and 
D. Belitz, Phys. Rev. Lett. 76, 2571 (1996). 

[8] D. Boyanovsky and J. L. Cardy, Phys. Rev. B 26, 154 (1982). 

[9] K.G. Wilson and M.E. Fisher, Phys. Rev. Lett. 28, 240 (1972). 
[10] S. Coleman and E. Weinberg, Phys. Rev. D 7, 1888 (1973). 
[11] A.B. Harris, J. Phys. C7, 1671 (1974). 
[12] N. Read and S. Sachdev, Nucl. Phys. B316, 609 (1989). 
[13] See J. Kogut, Rev. Mod. Phys. 55, 775 (1980) for a review. 

[14] H. Hamidian, S. Jaimungal, GW. Semenoff, P. Suranyi and L.C.R. Wijewardhana, 
Phys. Rev. D 53, 5886 (1996). 

[15] H. Hamidian, GW. Semenoff, P. Suranyi and L.C.R. Wijewardhana, Phys. Rev. Lett. 
74, 4976 (1995). 

[16] K.G. Wilson, Phys. Rev. D 14, 2911 (1974). 

[17] J. Zinn- Justin, Nucl. Phys. B367, 105 (1991); K. Karkkainen, R. Lacaxe, P. Lacock 
and B. Peterson, Nucl. Phys. B415, 781 (1994). 

[18] B. Rosenstein, H.L. Yu and A. Kovner, Phys. Lett. B314, 381 (1981). 

[19] S.F. Edwards in Polymer Networks, edited by A.J. Chompff and S. Newman (New York, 
Plenum Press, 1971); P.G. De Gennes, Phys. Lett. A38, 339 (1972); G. Grinstein et al. 
in above; V.J. Emery, Phys. Rev. B 11, 239 (1975); S.F. Edwards and PW. Anderson, 
J. Phys. F5, 965 (1976). 

[20] H. Yamagishi, Phys. Rev. D 23, 1880 (1981). 

[21] Lubensky in Ref.0 and A. Weinrib and B.I. Halperin, Phys. Rev. B 27, 413 (1983). 

Figure caption: Fig. 1: Projection of the RG trajectories in the six-dimensional coupling 
constant space onto the g\ — g<i plane for No = 2, Np = 13 and e — — 0.1. The dashed 
and solid curves correspond to theories with and without quenched disorder, respectively. 
The straight dashed line is defined by Eq. (18) and the lowest curve is a runaway trajectory. 
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